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Abstract 

This paper provides the basis for new methods of inference for max-stable 
processes £ on general spaces that admit a certain incremental representation, 
which, in important cases, has a much simpler structure than the max-stable 
process itself. A corresponding peaks-over-threshold approach will incorporate 
all single events that are extreme in some sense and will therefore rely on a 
substantially larger amount of data in comparison to estimation procedures 
based on block maxima. 

Conditioning a process n in the max-domain of attraction of £ on being 
extremal, several convergence results for the increments of n are proved. In 
a similar way, the shape functions of mixed moving maxima (M3) processes 
can be extracted from suitably conditioned single events n. Connecting the 
two approaches, transformation formulae for processes that admit both an 
incremental and an M3 representation are identified. 
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1. Introduction 

The joint extremal behavior at multiple locations of some random process 
{i](t) : t G T}, T an arbitrary index set, can be captured via its limiting max- 
stable process, assuming the latter exists and is non-trivial everywhere. Then, 
for independent copies rji of 77, i G N, the functions b n : T — > R, c n : T — > (0, oo) 
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can be chosen such that the convergence 

£(t) = lim c n (t)(ma<x.TH(t) - b n (tj), t £ T, (1) 

holds in the sense of finite-dimensional distributions. The process £ is said to 
be max-stable and n is in its max-domain of attraction (MDA). The theory of 
max-stable processes is mainly concerned with the dependence structure while 
the marginals are usually assumed to be known. Even for finite-dimensional 
max-stable distributions, the space of possible dependence structures is un- 
countably infinite-dimensional and parametric models are required to find a 
balance between flexibility and analytical tractability [7, 21]. 

A general construction principle for max-stable processes was provided by 
[6, 26]: Let XlieN ^(Ui,s z ) be a Poisson point process (PPP) on (0, oo) x S with 
intensity measure u~ 2 du ■ v(ds), where (5, 6) is an arbitrary measurable space 
and v a positive measure on S. Further, let / : SxT — > [0, oo) be a non-negative 
function with J s /(s, t)u(ds) = 1 for all t £ T. Then the process 

£(*)= max i), teT, (2) 

is max-stable and has standard Frechet margins with distribution function 
exp(— 1/x) for x > 0. In this paper, we restrict to two specific choices for / and 
(5, (3, u) and consider processes that admit one of the resulting representations. 
First, let {W(i) : t £ T} be a non-negative stochastic process with EW(t) = 1, 
t £ T, and W(to) = 1 a.s. for some point to £ T. The latter condition means 
that W(t) simply describes the multiplicative increment of W w.r.t. the location 
to- For (S, being the canonical probability space for the sample paths of 
W and with f(w,t) = w(t), w £ S, t £ T, we refer to 

Z(t)=m.BxUiWi(t), t£T, (3) 

as the incremental representation of £, where {VFjjigN are independent copies 
of W. Since T is an arbitrary index set, the above definition covers multivariate 
extreme value distributions, i.e. T = {t±, . . . , t^}, as well as max-stable random 
fields, i.e. T = R d . 

For the second specification, let {F(t) : t £ M. d } be a stochastic process with 
sample paths in the space C(lR rf ) of non-negative continuous functions, such 
that 

Ej Rd F(t)dt = l. (4) 

With Si = (Ti,Fi), i £ N, in S = R d x C(R d ), intensity measure v{dt x dg) = 
dtPp(dg) and f((t,g),s) = g(s — t), (t,g) £ S, we obtain the class of mixed 
moving maxima (M3) processes 



£(£) = maxJ7jFj(i-Tj), t £ M. d . 



(5) 
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These processes are max-stable and stationary on W d (see for instance [27]). 
The function F is called shape function of £ and can also be deterministic (e.g., 
in case of the Smith process). In Smith's "rainfall-storm" interpretation [26], 
Ui and Tj are the strength and center point of the ith storm, respectively, and 
UiFi{t — Ti) represents the corresponding amount of rainfall at location t. In 
this case, £(t) is the process of extremal precipitation. 

When i.i.d. realizations rjx, . . . ,rj n of n in the MDA of a max-stable process 
£ are observed, a classical approach for parametric inference on £ is based on 
generating (approximate) realizations of £ out of the data r]i,...,r] n via com- 
ponentwise block maxima and applying maximum likelihood (ML) estimation 
afterwards. A clear drawback of this method is that it ignores all information 
on large values that is contained in the order statistics below the within-block 
maximum. Further, ML estimation needs to evaluate the multivariate densities 
while for many max-stable models only the bivariate densities are known in 
closed form. Thus, composite likelihood approaches have been proposed [20, 5]. 
In univariate extreme-value theory, the second standard procedure estimates 
parameters by fitting a certain PPP to the peaks-over-thresholds (POT), i.e., 
to the empirical process of exceedances over a certain critical value [18, 10]. Also 
in the multivariate framework we can expect to profit from using all extremal 
data via generalized POT methods instead of aggregated data. In contrast 
to the ML approach, in this paper, we assume that £ admits one of the two 
representations (3) and (5) and we aim at extracting realizations of the processes 
W and F, respectively, from single extreme events. Here, the specification of a 
single extreme event will depend on the respective representation. 
In [12], this concept is applied to derive estimators for the class of Brown- 
Resnick processes [4, 17], which have the form (3) by construction. With a(n) 
being a sequence of positive numbers with lim n _^ 00 a(n) = oo, the convergence 
in distribution 



to,ti, . . . ,tk £ T, k £ N, is established for rj being in the MDA of a Brown- 
Resnick process and with W being the corresponding log-Gaussian random 
field. A similar approach exists in the theory of homogeneous discrete-time 
Markov chains. For instance, [25] and [9] investigate the behavior of a Markov 
chain {M{t) : t £ Z} conditional on the event that M(0) is large. The resulting 
extremal process is coined the tail chain and turns out to be Markovian again. 
In this paper, the convergence result (6) is generalized in different aspects. 
Arbitrary non-negative processes {W(i) : t £ T} with EW(t) = 1, t £ T, 
are considered, and convergence of the conditional increments of ry in the sense 
of finite-dimensional distributions as well as weak convergence in continuous 
function spaces is shown (Theorems 2.1 and 2.2). Moreover, in Section 3, 
similar results are established for M3 processes (5) by considering realizations 




(6) 
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of 77 around their (local) maxima. Since one and the same max-stable process £ 
might admit both representations (3) and (5) we provide formulae for switching 
between them in Section 4. Section 5 gives an exemplary outlook on how our 
results can be applied for statistical inference. 

2. Incremental representation 

Throughout this section, we suppose that {£(£) : t £ T}, where T is an 
arbitrary index set, is normalized to standard Frechet margins and admits a 
representation 

t(t) = maxU i V i (t), teT, (7) 

where XligN ^u t 1S a PPP on (0, 00) with intensity u~ 2 du, which we call Frechet 
point process in the following. The {Vi}j g N ar e independent copies of a non- 
negative stochastic process {V(t) : t £ T} with EV(t) = 1, t e T. Note 
that (7) is slightly less restrictive than the representation (3) in that we do not 
require that V(to) = 1 a.s. for some to e T. For any fixed to e T, we have 

e(t)=max^(lp i=0 ^ {1) (t) + lp I= iV; {2) (t)), teT, (8) 

where {Pj}j g N are i.i.d. Bernoulli variables with parameter p = P(V(to) = 0) 
and the and vf 2 are independent copies of the process {V(t) : teT}, 
conditioned on the events {V(to) > 0} an(1 {^(^0) = 0}, respectively. 

Note that for k G N, to, . . . , t). e T, the vector H = (£(io)> . . . follows 
a (k + l)-variate extreme- value distribution and its distribution function G can 
therefore be written as 

G(x)=exp(-/z([0,x] c )), x€R fc+1 , (9) 

where \i is a measure on£= [0, oo) fc+1 \ {0}, the so-called exponent measure 
of G [21, Prop. 5.8], and [0,x] c = E \ [0,x]. 

The following convergence result provides the theoretical foundation for sta- 
tistical inference based on the incremental process V. 

Theorem 2.1. Let {rj(t) : t e T} be non-negative and in the MDA of some 
max-stable process £ that admits a representation (7) and suppose that r\ is 
normalized such that (1) holds with c n (t) = 1/n and b n (t) = for n e N and 
teT. Let a(n) — » 00 as n — > 00. For k e N and to, . . . , t}~ € T we have the 
convergence in distribution on R fc+1 



>7(*o) v(tk) 



rj(t ) > a{n) A (z, AV^) , n -)• 00, 



ya(n) ' rj(t )' ' ' ' ' rj(t ) 

where the distribution of AV"' 1 ' is given by 

P(AVW e dz) = (l-p)P(AvW e dz)E(V W {t )\AV^ = z), z > 0. 



(10) 
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Here, AV"W denotes the vector of increments [ ytrrj^ > • • • > yalfo) J w ^ re spect 
to to, and Z is an independent Pareto variable. 

Remark 2.1. Note that any process r\ that satisfies the convergence in (1) for 
a process £ with standard Frechet margins can be normalized such that the 
norming functions in (1) become c n (t) = 1/n and b n (t) = 0, n G N, t G T [21, 
Prop. 5.10]. 

Proof of Theorem 2.1. For X = (r/(to), • • • ,v(tk))> which is in the MDA of 
the random vector S = (£(io)> • • • )£(ift))i it follows from [21, Prop. 5.17] that 

lim mP(X/meB)=fi(B), (11) 

m— >oo 

for all elements B of the Borel c-algebra B(E) of E bounded away from {0} with 
fi{dB) = 0, where \i is defined by (9). For sq > and s = (si, . . . , s^) G [0, oo) fc , 
we consider the sets A So = (s , oo) x [0, oo) fc , A = Ai and B s = {x £ [0, oo) fc+1 : 
(jcW, . . . , &(*)) < x(°)s} for s satisfying P(AVW G d[0, s]) = 0. Then 

{n(t ) > s a(n), (v(h)/r}(t ), n{t k ) / n{t Q )) < s} = {X/o(n) G -B s fl ^ So }, 

since -B s is invariant under multiplication, i.e., B s = cB s for any c > 0. Thus, 
we obtain 



p ( v(to) > s a(n), (n(ti)/r](t ), rj(t k ) /rj(t )) < s rj(t ) > a(n] 
a{n) P(X/a(n) £ B S C) A(~) A So ) 



a(n)P(X/a(n) G A) 

(n^oo), (12) 



/x(s B n ,4 n A* 



M (A) 

where the convergence follows from (11), as long as n{d{B s n j4 n ^4s )} = 0. 
Let 

£«(t) = maxL/fV i (1) (i), t G T, (13) 

where SieN^t/.w ^ s a Poisson point process with intensity (1 — p)w _2 dti and 
let /Z 1 ) be the exponent measure of the associated max-stable random vector 
(£ (1) (*o),---,£ (1) (*&))• Then the choice A = (1 , oo) x [0, oo) fc guarantees that 
//(• n A) = /i^'f- PI ^4). Comparing the construction of £W in (13) with the 
definition of the exponent measure, we see that fj,^ is the intensity measure of 
the Poisson point process Ylien ( W (13 V- (1) (to) u (1) v ll) (t k )) on Hence, 



/■OO 

= / (1 -p)u~ 2 P(uV {1) (t ) > l)du 
Jo 

POO P 

= {l-p) u- 2 P(^«(to)€dy)d 

J0 Jfu-i.oo) 



/•OO 

yP(V ( - 1 \t )Gdy) = (l-p)EV ( - 1 \t Q ) = l, (14) 
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where the last equality follows from E7W(to) = EV(to)/(l — p). Furthermore, 
for s > 1 and s G [0,oo) fc with P(AV< 1 ) G d[0,s]) = 0, 



M (B 8 ninAJ/((l-j>)M(i)) 

u- 2 P[uV W {t ) >s , (uV W {h),... lU V W (t k )) <suV W {t 

o) d-u 



= / / u~ 2 P(V (1) (i ) G dy AV* 1 ) < s) P(AV^ < s)du 

J 7 [sow -1 , do) ^ ' 

= / / ysg 1 -P(V (1) (io) € dy AV (1) = z) P(AV^ g dz) 

V[0,s] i[0,oo) V 7 

= s 1 / E(V (1) (io) AV^ = z) P(AV (1) G dz). (15) 

J[0,s] V ' 

Equation (15) shows that the convergence in (12) holds for all continuity points 
s G [0,oo) fc of the distribution function of AV' 1 '. Since so > 1 was arbitrary, 
this concludes the proof. 

Remark 2.2. 1. If V^'(to) is stochastically independent of the increments 
AVf 1 ), we simply have P(AV< 1 ) G dz) = P(AVW G dz). 

2. If p = P(V(£o) = 0) = 0, the exponent measure [i of any finite-dimensional 
vector H = (£(io)> • • • , ? ■ ■ ■ j^fe G T, G N, satisfies the condition 
/U ({0} x [0, oo) fe ) = 0, and following Proposition 2.1, the incremental 
representation of 3 according to (3) is given by H_= maxjgN Ui • (1, AVj) T , 
where AVj, i G N, are independent copies of AV = AV' 1 '. 

3. If £ admits a representation (3), we have P(AVW G dz) = P(AV G dz), 
which shows that (6) is indeed a special case of Theorem 2.1. 

Remark 2.3. In the above theorem, the sequence a(n) of thresholds is only 
assumed to converge to oo, as n — > oo, ensuring that {rj(to) > a(n)} becomes 
a rare event. For statistical applications a(n) should also be chosen such that 
the number of exceedances 

n 

N(n) = J2 1 {v l (to)>a(n)} 



i=i 



converges to oo almost surely, where (rji)^ is a sequence of independent 
copies of r\. By the Poisson limit theorem, this is equivalent to the additional 
assumption that lim ra _ il00 a(n)/n = 0, since in that case nP(r](to) > a(n)) = 
n/a{n) — > oo, as n — > oo. 

Remark 2.4. [12] consider Hiisler-Reiss distributions [15, 16] and obtain their 
limiting results by conditioning on certain extremal events A C E. They 
show that various choices of A are sensible in the Hiisler-Reiss case, leading 
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to different iimiting distributions of the increments of rj. In case £ is a Brown- 
Resnick process and A = (l,oo) x [0,oo) fc the assertions of Theorem 2.1 and 
[12, Thm. 3.3] coincide. 

Example 2.1. (Extremal Gaussian process [23].) A commonly used class of 
stationary yet non-ergodic max-stable processes on M. d is defined by 

£(f) = maxUiYi{t), teR d , (16) 

where X^<=n * s a Frechet point process, Yi(t) = max(0, %(€)), and the 

Yi are i.i.d. stationary, centered Gaussian processes with E(max(0, = 1 

for all t £ R d [23, 3]. Note that in general, a t £ R d s.t. Yi(t ) = 1 a.s. does not 
exist, i.e., the process admits representation (7) but not representation (3). In 
particular, for the extremal Gaussian process we have p = P(V(to) = 0) = 1/2 
and the distribution of the increments in (10) becomes 



P(AVW g dz) = ~ E \Y(t ) | (Y(ti)/Y(to), Y(t k )/Y(t )) = z, Y(t ) > 
>((Y(t 1 )/Y(to),...,Y{t k )/Y(t )) Gdz|y(t ) >0). 



2 
• P 



While the Hiisler-Reiss distribution is already given by the incremental rep- 
resentation (3), cf. [16], other distributions can be suitably rewritten, provided 
that the cumulative distribution function and hence the respective exponent 
measure /i is known. 

Proposition 2.1. Let H = (£(*o)> ••• be a max-stable process on T = 

{to, . . . ,tk} with standard Frechet margins and suppose that its exponent mea- 
sure // is concentrated on (0, oo) x [0, cc) k . Define a random vector W = 
(W^ 1 ), . . . , iy( fc )) via its cumulative distribution function 

P(W < s) = n(B B nA), sG[0,cx)) fc , (17) 

where A=(l,oo)x [0,oo) k and B s = {x G [0, oo) k+1 : (x^, . . .,1^) < x^s}. 
Then, H allows for an incremental representation (3) with Wj, i £ N, being 
independent copies o/W. 

Proof. First, we note that (17) indeed defines a valid cumulative distribution 
function. To this end, consider the measurable transformation 

T : (0, oo) x [0, oo) fc -> (0, oo) x [0, oo) fc , (x , ...,x k )^ ( x , 

Then, T(B s nA) = (l,oo) x [0,s] and the measure fi T (-) = ^(r _1 ((l, oo) x • )) 
is a probability measure on [0,oo) fc . Since 

fx(B s nA) = ^(T- 1 ((l,oo) x [0,s])) = ^ T ([0,s]), 
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the random vector W is well-defined and has law [i T . 

By definition of the exponent measure, we have E = maxjgpjXj, where II = 
SieN * s a on E with intensity measure \i. Then, the transformed point 
process Til = X]ieN^(x <0) x (1) /x <0) x w /x (0) ) nas intensity measure 

A((c, oo) x [0, s]) = ii (T- 1 ((c, oo) x [0, s])) 

= fi(B s n ((c, oo) x [0, oo) fc )) = c - V(-Bs n A) 

for any c > 0, s £ [0, oo) k , where we use the fact that fj,, as an exponent 
measure, has the homogeneity property c _1 /i(dx) = /z(d(cx)). Thus, TIT has 
the same intensity as Ylien Vi.Wjli where SieN^* * s a Frechet point process 
and Wj, i £ N, are i.i.d. vectors with law P(W < s) = fJ,(B a n A). Hence, we 
have 

Z^^T-^(X?\XV/X?\...,XP/XP 
= maxT^ 1 ((Ui, WA) = max^Wi, 
which completes the proof. 

Example 2.2. (Symmetric logistic distribution, cf. [14]-) F° r T = {to, . . . ,t k }, 
the symmetric logistic distribution is given by 



P(£(*o) < xo, ■ ■ • ,£(*k) < Xk) = exp 



Q I . — (J 



1/9 



(18) 



for xo, . . . , Xk > and g > 1. Hence, the density of the exponent measure is 

/ k \ V9-( fc +l) / fe \ fc 

/i(da;o, • • • , dxfc) = ^ x^ 9 (ig - 1) J JJ x^Mxo . . . dx k - 

\i=0 ) \%=\ / i=0 

Applying Proposition 2.1, the incremental process in the representation (3) 
is given by 



k 



1/9-1 



P(W(h) < Sl ,... W(t k ) < s k ) = 1 + J2 S 



g 



2.1. Continuous sample paths 

In this subsection, we provide an analog result to Theorem 2.1, in which 
convergence in the sense of finite-dimensional distributions is replaced by weak 
convergence on function spaces. In the following, for a Borel set U C M. d , we 
denote by C(U) and C + (U) the space of non-negative and strictly positive 
continuous functions on U, respectively, equipped with the topology of uniform 
convergence on compact sets. 
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Theorem 2.2. Let K be a compact subset of M. d and {rj(t) : t G K} be a 
process with positive and continuous sample paths in the MDA of a max-stable 
process {£,(t) : t G K} as in (3) in the sense of weak convergence on C(K). In 
particular, suppose that 

— max?7j(-) A £(•), n — > oo. 
n i=i 

Let W be the incremental process from (3) and Z a Pareto random variable, 
independent ofW . Then, for any sequence a(n) of real numbers with a(n) —> oo, 
we have the weak convergence on (0, oo) x C(K) 



a[n 



rj(t ) > a{n)) ^ (Z,W(-)), 



as n tends to oo. 



Remark 2.5. Analogously to [28, Thm. 5], weak convergence of a sequence of 
probability measures P n , n G N, to some probability measure P on C(M rf ) is 
equivalent to weak convergence of Pnrj 1 to Prj 1 on C([-j, j] d ) for all j > 1, 

where rj : C(M. d ) — > C([— j, j] ) denotes the restriction of a function to the cube 
[~j)j] d - Hence the assertion of Theorem 2.2 remains valid if the compact set 
K is replaced by R rf . 

Proof of Theorem 2.2. As the process £ is max-stable and i] G MDA(£), 
similarly to the case of multivariate max-stable distributions (cf. Theorem 2.1), 
we have that 

lim uP(i)/u G B) = n(B) (19) 

u — ^oo 

for any Borel set B C C(K) bounded away from 0^, i.e., inf{sup sgA - f(s) : f G 
B} > 0, and with /J,(dB) = [7, Cor. 9.3.2], where fi is the exponent measure 
of £, defined by 

P(£(s) < Xj, s G Kj, j = 1, . . . ,m) 



exp 



-/i 



N / G C{K) : sup s6 ^. f(s) > xj for some j G {1, . . . ,m}| 



(20) 



for xj > 0, Kj C K compact. Thus, // equals the intensity measure of the 
Poisson point process XlieN ^U t WA-)- F° r z > and D C C{K) Borel, we 
consider the sets 

A z = {/ G C(K) : /(t ) > z} 
B D = {fG C{K) : /(-)//(to) G ^} 

and A = A\. Note that -Bd is invariant w.r.t. multiplication by any positive 
constant. Then, as VF(to) = 1 a.s., we have ^(A z ) = J z u~ 2 du = z^ 1 and for 
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s > 1 and any Borel set D C C(K) with P(W G 3D) = 0, by (19), we get 

P U(t Q )/a(n) > s , v(-)/v(to) G £> >j(to) > o(n)} 
_ a(n) P{r?(-)/q(n) £ A Bo n Pi A] 
a(n)P{rj(-)/a(n) G A} 

r 2 P{uiy(') G £ D }du 



= s 1 P{w(-)eD}, 

which is the joint distribution of Z and W(-). 

Example 2.3. (Brown- Resnick processes, cf. [4, 17].) For T = R d , d > 1, let 
{y(t) : t G T} be a centered Gaussian process with stationary increments, 
continuous sample paths and Y(to) = for some to £ Note that by [1, 
Thm. 1.4.1] it is sufficient for the continuity of Y that there exist constants 
C, a, 5 > 0, such that 

E \Y(s) - Y(t)\ 2 < - „ ° „,„ 

for all s,t G R d with < 5. Further let j(t) = E(Y(t)-Y(0)) 2 and a 2 (t) = 

~E(Y(t)) 2 , t G denote the variogram and the variance of Y, respectively. 
Then, with a Frechet point process XlieN ^u, and independent copies Yi of Y, 
i G N, the process 

f(t) = max£7 i exp(F i (i) -a 2 (t)/2) , t G M d , (21) 

is stationary and its distribution only depends on the variogram 7. Comparing 
(21) with the incremental representation (3), the distribution of the increments 
is given by the log-Gaussian random field W(t) = exp (Y(t) — <j 2 (t)/2) , t G M. d , 
and Theorem 2.2 applies. 

3. Mixed moving maxima representation 

A large and commonly used class of max-stable processes is the class of M3 
processes (5). Let 

110 = 'Yl S (Ui-Ti,F i ) ( 22 ) 

be the corresponding PPP on (0, oo)xM d xC(M d ) with intensity u~ 2 dudt P_p(d/). 
In the sequel, M3 processes are denoted by 

M(t) =maxU i F i (t-T i ), t G M. d . 
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The marginal distributions of M are given by 



P(M(t ) < s ,...,M(t k ) < s k ) 



n ( { (u,t,f) : maxu/(i, - t)/si >U I - 



exp 



f [ nu«(/(t,-t)/ 8 ,)d*P F (d/) ) , 



, tfc G R , s ,..., a* > 0, k G N. 



(23) 



In Section 2, we were interested in recovering the incremental process W from 
processes in the MDA of a max-stable process with incremental representation. 
In case of M3 processes, the object of interest is clearly the distribution of the 
shape function F. Thus, in what follows, we provide the corresponding con- 
vergence results for processes r\ in the MDA of an M3 process. We distinguish 
between processes on R d with continuous sample paths and processes on a grid 
(Z d ). The main idea is to consider n in the neighborhood of its own (local) 
maximum, conditional on this maximum being large. 



3.1. Continuous Case 

Let {r](t) : t G R rf } be strictly positive and in the MDA of a mixed moving 
maxima process M in the sense of weak convergence in C(R rf ). We assume that 
7] is normalized such that the norming functions in (1) are given by c n (t) = 1/n 
and b n (t) = 0, for any n £ N and t G K rf . Further suppose that the shape 
function F of M is sample-continuous and satisfies 



F(0) = A a.s., 
F(t) G [0, A) Vt G 



a.s. 



for some A > and 



l E {s F(t °-4 dt<0 ° 



(24) 



(25) 



for any compact set K C Mr. Under these assumptions, there is an analog 
result to Theorem 2.2. 

Theorem 3.1. Let Q,K cR d be compact such that dQ is a Lebesgue null set 
and let 

t q : C(Q) ^R d , f 1 ^ inf ( arg max f(t) ) , 

\ teQ J 

where "inf" is understood in the lexicographic sense. Then, under the above 
assumptions, for any Borel set B C C{K) with P(F/X G dB) = 0, and any 
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sequence a(n) with a(n) — > oo as n — >■ oo, we have 
lim limsupP {^(tq^q) + •) /v(.tq{v\q)) G 5 

max7?(i) = max n(t), maxr?(i) > a(n) > = P { F(-)/X G £>), 



{ojeL/-: 

compact 



where © denotes morphological dilation. 

The same result holds true if we replace limsup,^^ by liminfn^oo. 

Proof. First, we consider a fixed compact set L C M. d large enough such that 
K U {0} C L and define 



A L = {f £ C(Q®L) : max/0) > 1, max/0) = m &x /(t) 
and 

C B = {/GC(Q©L): /(tq(/| )+ -)//(ro(/|o))€B} 

for any Borel set B C C(K). Note that C# is invariant w.r.t. multiplication by 
any positive constant. Thus, we get 

P{v{tq(v\q) + ■) / V(tq(v\q)) G B max7?(t) = max rj(t) > a(n)\ 

= P (77/0(71) G C B I r//o(n) G A L } 
_ a(n) P {77/0(71) £ C B , 77/0(71) £ A L ) 
a(n) P {r?/a(n) G A L ) 

By [8, Cor. 9.3.2] and [21, Prop. 3.12] we have 

limsupitP(?7/u eC)< m(C)j C C C(Q © L) closed, 

u— >oo 

liminf uP(??/u eO)> ^t(O), O C C(Q © L) open, 

u— ¥00 

where C and O are bounded away from 0^. Here, [i is the intensity measure of 
the PPP Ylien $UiFA ■ -TA restricted to C{Q@L). Thus, by adding or removing 
the boundary, we see that all the limit points of Equation (26) lie in the interval 



(26) 



ii{C B n Al) - Kd(C B n A L )) n(C B nA L )+ii(d(C B nA L )) 
H{A L ) + n(8A L ) ' n(A L ) - n(8A L ) 

We note that Al is closed and the set 
A* L = ^f eC{Q@L): 

TQ{f\Q)€Q°, max/(t) >max{l,/(t)} Vi G Q © L \ {tq(/|q)} 



(27) 
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is in the interior of Al (Lemma A.l). Hence, we can assess 
KdA L ) < /x({/ G C(Q © L) : max/(i) = 1}) 



+ M {/GC(Q0L): tq(/|q) G <9Q} 



U < / G C(Q © L) : argmax/(t) is not unique 
[ teQ®L " J 

n <j / G C(Q © L) : max/(i) = max f(t) > 1 

oo 

< + / / u~ 2 du dt 



+ 



1) 



9Q 



/ / u~ 2 P<jumaxF(i -x) > lldiidx. (28) 



Here, the equality /i({/ G C(Q © L) : max ie q/(t) = 1}) = holds as 
max tG Q M(t) is Frechet distributed (cf. [8, Lemma 9.3.4]). Since dQ is a 
Lebesgue null set, the second term on the right-hand side of (28) also vanishes. 
Thus, 

K 9A l) < / / u ~ 2 p \umaxF(t - x) > 1 1 dudx =: c(L). (29) 

Now, let 5 C C(if) a be Borel set such that P(F/A G dB) = 0. For the set 
Cbi we obtain that the set 

C* B = if e C(Q © L) : argmax/(t) is unique, G £° 

I /eQ J{ t QU\q)) 



C* B U < / G C(Q © L) : argmax/(t) is not unique 
{ t<EQ J 



is in the interior of Cb and that the closure of Cb is a subset of 

■) 

U {/ G C(Q © L) : /(r Q (/| Q ) + -)//M/|q)) G 55} 
(Lemma A. 2 and Lemma A. 3). Thus, by (29), we can assess 

n(d(c B n a l )) < n(dA L ) + n(dc B n a l ) 

< c(L) + / / u~ 2 P jnmaxF(t - x) > 1 1 dudx 

Jri\(QS)L) Jx- 1 { *oGQ J 

f oo 



+ / / u~ 2 P(F/A G dB)dudt = 2c(L). (30) 
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Furthermore, we get 
fi(C B n A L ) 
= j j u- 2 p{f(-)/\£ B^dudto 

/ u~ 2 P < u max F(tn — x) > 1, 
Jx-i { *o€Q 

F ((tq(F(- - x)\ Q )) + ■ - x) I maxF(i -x)eB, 



+ 



d \(Q(BL) 



t &Q 

F(t - x)/maxF(t - x) < 1 Vt G Q®L \dudx. (31) 

The second term in (31) is positive and can be bounded from above by c(L). 
Setting B = C(K), ^(Al) can be expressed in an analogous way. Now, we plug 
in the results of (29), (30) and (31) into (27) to obtain that all the limit points 
of (26) are in the interval 

A • |Q| • P {F(-)/X £ B}-2c(L) A - \Q\ -P{F(-)/X £ B) + 3c(L) 
A • \Q\ + 2c(L) ' X-\Q\-c(L) 

Finally, we note that c{L) can be bounded from above by 

E < maxF(to — x) > dx, 

d \(Q®L) L *o6Q J 

which vanishes for L /* M. d because of assumption (25). This yields the assertion 
of the theorem. 

We conclude the treatment of the continuous case with an example of a 
process r/ that allows for an application of Theorem 3.1. As ry will be composed 
of a (locally) finite number of shape functions from the M3 construction in (5), 
r\ may directly model rainfall data and has therefore the potential for various 
practical applications. 

Example 3.1. Let {F(t) : t G M. d } be a random shape function as defined in 
(4). For c,e > let n c , e = E iS N be a PPP on (0,oo) x R d x C(R d ) 
with intensity 

cl{„> e }U'" 2 d'udt P F {df)- 
and, for k > 0, define a process M = M c ,e,n by 

M(-) = k V max uf(- -t). 



Then, the following statements hold. 
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1. M is in the MDA of the M3 process M associated to F in the sense of 
finite-dimensional distributions. 

2. If F satisfies (25), then M is in the MDA of M in the sense of weak 
convergence on C(IR d ). 

For a proof of this example, the reader is referred to Appendix B. 

3.2. Discrete Case 

Theorem 3.1 allows for estimation of F if the complete sample paths of 
7] are known, at least on a large set Q © L C For many applications, 

this assumption might be too restrictive. Therefore, we seek after a weaker 
assumption that only requires to know r/ on a grid. This needs a modification 
of the underlying model leading to a discretized mixed moving maxima process. 

Let {F(t) : t G be a measurable stochastic process with values in [0, oo) 
and 

J>F(*) = 1. (32) 

Further, let n 0j discr = Ylien &(Ui.Ti,FA be a Poisson point process on (0, oo) x 
7L d x [0,oo) z<i with intensity u~ 2 du 5%d (di) P^(d/). Then, the discrete mixed 
moving maxima process M^ SCT is defined by 

M discr (i) = max UiFi{t -7}), t G Z d . (33) 

The process Mdi scr is max-stable and stationary on 7L d and has standard Frechet 
margins. 

Let {rj(t), t G Z d } be in the MDA of a discrete mixed moving maxima 
process Mdi SC r in the sense of convergence of finite-dimensional distributions 
with norming functions c n (t) = l/n and b n (t) = in (1), n G N and t G Z d . 
Furthermore, we assume that the shape function F satisfies (24) with M. d being 
replaced by Then, analogously to Theorem 3.1, the following convergence 
result can be shown. 

Theorem 3.2. Under the above assumptions, for any k G N, k + 1 distinct 
points to, . . . , tk G Z d , any Borel sets B%, . . . , C [0, oo) such that 

P (OF(ii)M, • • • , F(t k )/X) G 8(B 1 x---xB k )}=0, 

and any sequence a{n) with a(n) — > oo as n — )• oo ; i£ holds 

lim lim P {77(^0 + *i)M*o) S B„ i = 1, . . . , k I 

{0}eL »Woo 

compact v(to) = maWto + 1), r?(t ) > o(n)} 



= P{W£B;, i = l,...,k}. 
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4. Switching between the different representations 

In the previous sections we analyzed processes that admit the incremental 
representations (3) or (7) and, on the other hand, processes of M3 type as in (5). 
We show that under certain assumptions, we can switch from one representation 
to the other. 

4.1. Incremental representation of mixed moving maxima processes 

We distinguish between M3 processes with strictly positive shape functions, 
for which we can find an incremental representation (3), and general non- 
negative shape functions, for which only the weaker representation (7) can 
be obtained. 

4.1.1. Mixed moving maxima processes with positive shape functions 

Theorem 4.1. Let M be an M3 process on M. d as in (5) with a shape function 
F with F(t) > for all t e R d . Then M admits a representation (3) with 
to = and incremental process W given by 

P(W€L)= [ [ l m ._ t)m _ t)eL} f(-t) dt P F (df), L G B(C+(R d )). 

JC+(R d ) JR d 

(34) 

Proof We consider the two Poisson point processes on (0,oo) x C + (M. d ) 

n i = fi(UjFi(—Ti),Fi( — Ti)/Fi(—Ti)) i (35) 

transformation of IIo in (22), and 

n 2 = z2 S (Ul,Wi(-)), ( 36 ) 

with Wi, ieN, being independent copies of W, and with Y^ien^U'- being a 
Frechet point process. Then the intensity measures of III and II2 satisfy 



EIIi([z,oo) x L) 



C+(R d ) 



U 



-2 



l{uf(-t)>z}l{f(-t)/f(-t)eL} dudt Pj?(d/) 



= z- 1 I [ l {/( ._ t)//( _ t)er} /(-t) dt P F (df) 

JC+(R d ) JR d 

= z' x Y{W e L) 
= En 2 ([z,oo) x L), 

L e B(C + (M> d )), z > 0, and hence IT = II2. The assertion follows from the fact 
that M is uniquely determined by 111 via the relation M(t) = max^ v g ^ Ui vg(t), 
t e R d . 
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While the definition of W in (34) is rather implicit, in the following, we 
provide an explicit construction of the incremental process W, which can also 
be used for simulation. To this end, let XlieN ^f" be a Frechet point process 
and let the distribution of (5, G) G C+(R d ) x R d be given by 

P((S,G)e(BxL)) (37) 

= / / l s eBl feL YW^r ds ( [ f(r)dr) P F (df) 
Jc+(R") Jr" J f{r)dr \J ) 

= [ [ l seB l feL f(-s)dsP F (df), 

JC+(R d ) JM d 

B G B d , L G B(C + (R d )). In other words, P G (d/) = (//(r)dr) P F (df) and, 
conditional on {G = /}, the density function of the shift S is proportional to 
/(— )• Putting W(-) = G(—S)/G(-S), equation (34) is satisfied and with i.i.d. 
copies Wi, i G N, of W, we get that maxjgN U"Wi(-) is indeed an incremental 
representation (3) of the mixed moving maxima process M. 

Remark 4.1. [MS representation of Brown-Resnick processes, cf. [17].) We 
consider the following two special cases of mixed moving maxima processes: 

1. Let £ £ M. dxd be a positive definite matrix and let the shape function 
be given by F(t) = (2ir)- d / 2 \T l \- 1 / 2 exp {-^Z^t}, t G R d . Then, M 
becomes the well-known Smith process. At the same time, by (37), S ~ 
iV(0, S) and G = F. Thus 

Y{t) = exp {-i(t - 5) T S- 1 (t - S) + i^S- 1 ^} 
= exp {-i^E-^ + i T E-^} . 

Since E(t T S _1 5) 2 = t T S _1 t, M is equivalent to the Brown-Resnick pro- 
cess in (21) with variogram j(h) = h T Y,~ 1 h. 

2. For the one-dimensional Brown-Resnick process £ in (21) with variogram 
j(h) = \h\, i.e., Y is the exponential of a standard Brownian motion with 
drift — \t\/2, [11] recently showed that the M3 representation is given by 
{F(t) : t £ 1} = {Y(t) | Y(s) < Vs G R : t G R}, i.e., the shape 
function is the exponential of a conditionally negative drifted Brownian 
motion. Having these two representations, it follows that the law of the 
conditional Brownian motion F, re-weighted by J F(t)dt and randomly 
shifted with density F{—-)/ J F(t)dt, coincides with the law of Y. 

44.2. Mixed moving maxima processes with finitely supported shape functions 

Let M be an M3 process on R d as in (5). In contrast to Section 4.1.1, where 
the shape functions are required to take positive values, here, we allow for 
arbitrary shape functions with values in [0, oo). 
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Theorem 4.2. The M3 process M as in (5) allows for an incremental repre- 
sentation of the form (7), with incremental processes V$ given by 

Vi(-) = Fi(- - Ri)/g(Ri). 

Here Ri, iGN, are i. i. d. copies of a random vector R with arbitrary density g 
satisfying g(t) > for all t G M. d , and F{, i G N, are i.i.d. copies of the random 
shape function F. 

Proof. With YlieN fat being a Frechet point process, we consider the process 
M(t) = maxUiFi(t- Ri)/g(Ri), t£R d , 

iGN 

which clearly is of the form (7). Then, 
P(M(«o) < s ,...,M(t k ) < s k ) 

= expf-/ [ max(/(^-t)/( 5 (t) Sz ))ff(t)dt P F (d/) ] 

\ JC(R d ) JR d 1=0 J 

= exp [ - ! ( max(/(t, - t)/ Sl )) dt P F (d/) j . 
\ JC(R d ) JR d 1=0 J 

The right-hand side coincides with the marginal distribution of M, which is 
given by (23). This concludes the proof. 

Decomposing V as in (8) with to = 0, we obtain the equality in distribution 
F (D(.) R _ R)/g{R) | _ R G S upp(F)). 
Applying Theorem 2.1 yields 
P(AV (1) G dz) 

/•oo 

= P(F(-R)/g(R) > 0) • / yP (V {1) (0) G dy, AV (1) G dz) 
g(s)ds P F (df) 



d ) J-B\ipp(J) 

oo r r 

V I / 1 /(-t)/9(t)Gdyl(/(t ! ~t)//(-t))f =1 edz 

iC(E d ) 7-supp(/) 

(/-suppm^d^'dt P F (d/)dy 



</(-*) da P F (d/) 

d ) ^supp(/) 



d ) Jsupp(f) 



/(*) 1 (/(t«+t)//(t))f =I ed. (/supp(/) 5(-«)ds) dt P F (d/). 

(38) 
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If the shape function F is deterministic, the right-hand side of (38) simplifies to 
/supp(/) /( t ) 1 (/(t ! +t)//(t))tiedzdt, i.e., the asymptotic conditional increments of 
7] G MDA(M) can be seen as a convolution of the shape function's increments 
with a random shift, whose density is given by the shape function itself. Note 
in particular, that this distribution is independent of the choice of the density 
g in Theorem 4.2. 

Remark 4.2. Section 4.1.1 considers the subclass of M3 processes with strictly 
positive shape functions and provides an incremental representation as in (3), 
which is nicely related to the conditional increments of r\ due to the property 
W(0) = 1. Section 4.1.2 applies to arbitrary M3 processes but only yields an 
incremental representation as in (7), for which the incremental process V does 
not directly represent the conditional increments of n. 

4.2. Mixed moving maxima representation of the incremental construction 

Theorem 4.3. Let Ylien^ be a Frechet point process and let Wi, i £ N, be 
independent copies of a non-negative, sample- continuous process {W(t), t € 
M. d }, satisfying 

lim W{t) = a.s., 

\\t\\— Kx> 

E W{t) = 1 for all t € R d , 

and ~E {m&xteK W(t)} < oo for any compact set K C M. d . 

Furthermore, let W be Brown-Resnick stationary, i.e., the process £, defined by 

Z(t) = maxUiWi(t), teR d , 

is stationary with standard Frechet margins. Then, the following assertions 
hold: 

1. The random variables 

Ti = inf < arg sup Wi (t) > and ji = sup Wi (t) 
I teR d J teR d 

are well-defined. Furthermore, X^gn &(Ui~fi,T l ,w l (-+T l )/'y i ) * s a Poisson point 
process on (0, oo) x x C(IR d ) with intensity measure ^>(du, dt,df) = 
cu~ 2 du dt P p(df) for some c > and some probability measure Pp. 

2. £ has an M3 representation with Pj?(d/) = P^.(cd/) being the probability 
measure of the shape function F. The constant c > is given by 



/(t)P # (d/)dt (39) 

C(R d ) ' I 
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and the probability measure P s is defined by 

p , A s _ J °°yP(W(- + T)/y £ A, r £ K \ j = y) P 7 (dy) 
F[ ' ^yP(r£K\ 1 = y)P 1 {dy) 



for any Borel set A C C(K ) and any compact set K C Mr, where r and 
7 are defined as 73 and ji, respectively, replacing W\ by W, and P 7 is the 
probability measure belonging to 7. 

Proof. 1. Analogously to the proof of [17, Thm. 14]. 

2. From the first part it follows that 

^0 = ^ , %,7i/c,T.,oH / .(-+T.)/7.) 

is a PPP with intensity measure ^>Q(du, dt, df) = u~ 2 du x dt x P^?(d/) 
where Pi?(d/) = P^,(cd/). Hence, $0 is of the same type as Ilo from the 
beginning of Section 3 and 

£(*)= max iGM d , 
(j/,s,/)e*a 

is a mixed moving maxima representation. The integrability condition (4) 
follows from the fact that £ has standard Frechet marginals. Thus, 



! ( cf(t)Pp(df)dt = l, 

JW JC(R d ) 



') 

which implies (39). In order to calculate P^,, let A £ B(C(M d )) and 
K £ B d be compact. The first part of this Theorem implies that 

*([1, 00) x K x A) = c ■ \K\ - Pp(A). 

Therefore, 

p - (A) = *([l.°°)*gxA) m 

and both the enumerator and the denominator are finite. For the enumer- 
ator, we get 

*([l,oo) x K x A) 

POO P OO 

= / u- 2 / P(W(- +r)/ 7 G A, r G K I 7 = y) P 7 (dy)dn 

JO Ju- 1 
roo roo 

= / u~ 2 du-P(W(- + T)/y £ A, T eK\j = y) P 7 (dy) 

/■oo 

= / y P(W(- + T)/y £A, r£K\ 1 = y) P 7 (dy). 
Jo 
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Thus, by (40) 



J™yP(W(- + T)/yeA, T€K\ 1 = y) P 7 (dy) 
J™yP(TGK\ 7 = y) P 7 (dy) 



which completes the proof. 



5. Outlook: Statistical applications 



In univariate extreme value theory, a standard method for estimating the 
extreme value parameters fits all data exceeding a high threshold to a certain 
Poisson point process. This peaks-over-threshold approach has been generalized 
in [22] to the multivariate setting. Therein, generalized multivariate Pareto 
distributions are obtained as the max-limit of some multivariate random vector 
in the MDA of an extreme value distribution by conditioning on the event that 
at least one of the components is large. Conditioning on the same extremal 
events, the recent contribution [13] analyzes the asymptotic distribution of 
exceedance counts of stationary sequences. 

Here, we have suggested conditioning a stochastic process rj(t) : t G T} in 
the MDA of a max-stable process {£(t) : t £ T} such that it converges to 
the incremental processes W in (3) or the shape functions F in (5). In this 
section we provide several examples how these theoretical results can be used 
for statistical inference. The approach is based on a multivariate peaks-over- 
threshold method for max-stable processes, though the definition of extreme 
events differs from that in [22, 13]. 

In the sequel, suppose that rji, . . . ,rj n , n 6 N, are independent observations of 
the random process rj, already normalized to standard Pareto margins. 

5.1. Incremental representation 

For a max-stable process £ that admits an incremental representation 



as in (3), the statistical merit of the convergence results in Theorem 2.1 and 
Theorem 2.2 is the "deconvolution" of U and W which allows to substitute 
estimation of £ by estimation of the process W. As only the single extreme 
events converge to W, we define the index set of extremal observations as 



for some fixed to £ T. The set {??i(-)/ r /i(*o) : * £ h(n)} then represents 
a collection of independent random variables that approximately follow the 
distribution of W. Thus, once the representation in (41) is known, both 
parametric and non-parametric estimation for the process W is feasible. For 
statistical inference it is necessary that the number of extremal observations 
|Ii(n)| converges to oo, as n — > oo. This is achieved by choosing the sequence 
of thresholds a(n) according to Remark 2.3. 



£(t) = maxUiWi(t), t£T, 



(41) 



h(n) = {i € {1, . . . n} : 7^ (to) > a(n)} 
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Example 5.1. (Symmetric logistic distribution, cf. Example 2.2.) The de- 
pendence parameter q > 1 of the symmetric logistic distribution (18) can 
be estimated by perceiving the conditional increments of rj in the MDA as 
realizations of W and maximizing the likelihood 

P(W(t 1 )€dai,...W(t k )€da k \q) 

= (i + Eti s~ q ) 1AMfe+1) (uU^ - 1)) nto ^ ^. 

Example 5.2. (Brown- Resnick processes.) Recall that the Brown- Resnick pro- 
cesses in Example 2.3 admit a representation (3) with log-Gaussian incremental 
process W(t) = exp {Y(t) — <7 2 (i)/2| , t G R d . Hence, standard estimation pro- 
cedures for Gaussian vectors or processes can be applied for statistical inference. 
[12] explicitly construct several new estimators of the variogram 7 based on the 
incremental representation, which also covers Hiisler-Reiss distributions, and 
they provide some basic performance analyses. 

5.2. Mixed moving maxima representation 

Similarly, in case of the mixed moving maxima representation 

M(t) = maxUiFi(t -Si), t G R d , (42) 

the convergence results of Theorem 3.1 can be used to estimate F (or Fx = F/\) 
on some compact domain K instead of estimating M directly. Here, the index 
set T of the observed processes {rji(t) : t G T}, i = 1, . . . , n, can be identified 
with Q®L from Theorem 3.1. The set L should be sufficiently large such that it 
is reasonable to assume that the components {UiFi(- — Si) : Si ^ Q© L] hardly 
affect the process M on Q®K (that is, \i{C B n A L )/fi(A L ) m P(F(-) - A G B) 
in the proof of Theorem 3.1). At the same time, a large set Q leads to a rich 

set of usable observations F^' = ^(tq^) + O/^^qC 7 /*))) * e ^(n), where 

I2M = U 6 {1, . . . , n} : maxr/j(t) = max rji(t) > a(n) > . 
{ teQ t£Q®L ' J 

The resulting processes F[ , i G hin), can be interpreted as independent 
samples from an approximation to F\. This approach can be expected to be 
particularly promising in case of F having a simple distribution or even being 
deterministic. 

Example 5.3. (MS processes with deterministic shape functions.) Some ex- 
amples of mixed moving maxima processes have already been analyzed for 
statistical inference by [8] who use normal, exponential and t densities as shape 
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functions. More precisely, they consider M3 models with 

F 1 (t)=exp(-^1, A = -g=, (43) 



2 J V2tt 
iq(i) = exp{-/3|i|}, A = |, (44) 

and + ">°' (45) 

all parametrized by f3 > 0. [8] introduce consistent and asymptotically normal 
estimators based on the interpretation of f3 as a dependence parameter. From 
the samples F^\ i G l2( n ), we get a new estimator 



Lvf(') 



\Hn„ : _ t 

for Applying this estimator, f3 can be estimated by a least squares fit of 
(43)-(45) to F\ at some locations ti,...t m £ K. Note that in case of the 
normal model (43) and the exponential model (44), the logarithm of the shape 
function F\ depends linearly on /3 2 and /3, respectively, and log F\ can be fitted 
by ordinary least squares. 

Example 5.4. (Brown- Resnick processes.) The mixed moving maxima rep- 
resentation can also be employed for estimation of Brown-Resnick processes 
although the distribution of F is much more sophisticated than the one of W 
in the incremental representation (cf. [11, 19]). A relation between the shape 
function F and the variogram 7 of the Brown-Resnick process can be obtained 
via the extremal coefficient function &(■)■ For a stationary, max-stable process 
£ with identically distributed marginals, [24] defined the extremal coefficient 
function 6 via the relation 

P(C(0) < u, £(/i) < u) = P(£(0) < u) e( - h \ h G R d . 

For mixed moving maxima processes, we have 

m = E / m) V Fit + H)} * = E ^'f V /}^" )}dt (46) 

jR d & J R d ri{t) at 

and, at the same time, for Brown-Resnick processes [17], 



6(h) = 2$ (VtW/2) , 



(47) 
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where <& is the standard Gaussian distribution function. Identifying (46) with 
(47) and plugging in the samples F^ , i G 1% (n) , we get the variogram estimator 

2 



2$ _1 



V / F^(t)VF^>(t + h)dt 



\ 



2 £ /^i (*) d * 
'if 



where K is a large set such that K,K + h C if 



\ 1 



/ J 



Appendix A. Auxiliary Results for the Proof of Theorem 3.1 

Lemma A.l. Al is closed. The set A* L is in the interior of A^. 

Proof. The first assertion is obvious. For the second one, let /* G A* L . 
Then, we have /*(tq(/*|q)) =: a > 1. Furthermore, there is 5 > such that 
B S (tq(J*\q)) = {t G M d : \\t - tq(J\q)\\ < 5} G Q° and we have 



P := sup f*(t) - max/*(i) < 0. 

teQ®L\B s (T Q (f*\ Q )) *eQ 



(48) 



Now, we choose s < min{^-, ^} and show that B e (/*) = {/ G C(Q L) : 



11/ - /*||oo < e} C ^l- This holds, as for any / G B e (f*), we have 



/(r Q (/| Q ))>/(TQ(r| Q ))>a-e> 



1 + a 



> 1 



and max/(i) < max /(t)=max< max f(t), max/(i) 



te<2 



t€Q(BL\Q° 



teQ 



< max < /3 + a + e, max/(i) > < max < a — e, max/(i) 



max f(t), 



which means equality. 

Lemma A. 2. The set C* B is in the interior of Cb- 

Proof. Let /* G C B . Then, t* = arg max t€ Q f*(t) is well-defined and neces- 
sarily, as / > 0, 

a:=f*(t*)e(0,\\r\U. (49) 
Since f*(t* + ■)//*(**) G tnere is some e > such that 



{/ G C(K) : ||/*(t* + •)/ /*(<*) " /|| < 4 



< e \ C 5. 



(50) 



Furthermore, /* is uniformly continuous on the compact set Q © L, i.e. there 
exists some 5 > such that 



sup 

s.teQffiL, ||s-t||<5 



irw-rcoi<o«- 



(51) 
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Then, as argmax tg g /*(£) is unique, we have that 



p:= max f*(t) -/*(**) G [-a, 0). 



(52) 



Chooser* < min{^},fp^-}. We will show that B e * a (f*) = {/ G C(QffiL) : 
11/ — /*l|oo < £*a} C Cb- To this end, let /o G B e * a {f*). Then, because of 
Equation (52) and e*a < we have that \\to — t*\\ < 5 for to = tq(/o|q)- 
Therefore, 



sup 

t&K 

< sup 

t&K 



/*(** + •) /o(*o + 



/*(**) 

/*(t* + ■ 



/o(<o) 



+ sup 

teK 



f*(t*) 

F(to + 



/*(**) 
/o(*o + 



+ sup 



r(*o+o r(to + 



/*(**) 



/o(*o 



/o(*o) 



< 



/o(*o) 
e* ||/* 



3 + 1 - e* 



+ 



a 



(53) 



where we used Equation (51) and the fact that /o G B e * a (f*). Equations (49) 
and (52) and the choice of e* yield that e*/(l-e*) < (e*||/*|| 00 )/((l -e*)a) < 
2e* 1 1/*| loo/a < e/3, i.e. each summand on the right-hand side of (53) is smaller 
than e/3. Thus, f (t + -)//o(*o) G {/ G C{K) : ||/*(t* + •)//*(**) - /IU < 
e} C B by Equation (50) and /o G Cb. 

Lemma A. 3. T/ie closure of Cb is a subset of 



B* U < / G C(Q © L) : argmax/(i) is not unique > 
[ " *eQ J 

U {/ G C(Q © L) : /(r Q (/| Q ) + •) /f(r Q (f\ Q )) G Si?} . 

Proof. Let {/ n } C C# be a sequence converging uniformly to some /* G 
C7(Q © L). We have to verify that F(tq(J*\q) + •) / 'f '(tq(/*|q)) e B U OB 
if argmaxjgQ /*(t) is unique. Analogously to the proof of Lemma A. 2 we can 
show that for any £2 > there is some e\ > such that 



11/ - /*||oo,QeL < £1 
/(tq(/|q) + /*(tq(/*|q) + -) 



/(tq(/|q)) 



/*M/*|q)) 



< £ 2 - 



oo,_ft" 



Thus, /n(r Q (/n|o) + -)//«(ro(/nlo)) converges to /*(t q (/1 q )+-)//*(tq(/*|q)) 
in C(A). Hence, as BUdB is closed, /*(r Q (/*| Q ) + •)//* (tq(/*|q)) G #Ud£. 
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Appendix B. Proof of Example 3.1 

Proof of Example 3.1. Let Mj, j G N, be independent copies of the process 
M and consider 

1 n ~ 
M n (•) = — maxMi(-). 
cn t=l 

Further, suppose that L C M rf is an arbitrary compact set. Note that by Remark 
2.5 it suffices to show weak convergence of M n — > M, n — > oo, on C(L). 

To prove the first assertion, note that, for to,. . . ,tk 6 sq, . . . ,Sk > 0, 
k 6 N, we have 



P(M n (t ) < s ,...,M n (t fc ) < s k 

U<minf =0 cn Sl ' P j N («,(,, j j . uiaj 



l«<mintn ens, • P 1 n c,e ( \ (u,t, /) : maxu/(i/ - t)/(cnsi) >lM = 



= l K <minf_ cn Si - exp I -n f f min ( - , max — — - j c dt P F (d/) J 
— > exp ( - / / max(/fa - t)/a { ) dt P F (d/) j , (54) 

\ JC(IR d ) JR d «=° / 

as n — )• oo, where the convergence holds due to monotone convergence. The 
right-hand side of (54) coincides with the marginal distribution of M (cf. (23)). 

For convergence of M n to M in the sense of weak convergence in C + (L) 
endowed with the topology of uniform convergence, it remains to show that the 
sequence of restricted processes {M n \i : n 6 N} is tight. To this end, by [2, 
Thm. 7.3], it suffices to verify that for any e > 0, rj 6 (0, 1), there exist 5 > 0, 
no G N such that 

pi sup \M n (s) - M n (t)\ >e\ <rj, n > n . 
[\\s-t\\<6 J 

By Equation (25), we can choose R > such that 

f E[supF(t-s) J ds < ^, (55) 

where Br(0) = {x £ R d : \\x\\ < R}. Furthermore, (25) implies that 
E (sup te ^ F(s)) < oo for any compact set K C Therefore, as each realiza- 
tion of F is uniformly continuous on B R+d ^(0), where d(L) = sup Si S2eL \ \s\ — 
S2II denotes the diameter of L, dominated convergence yields 

limE I sup \F(s) - F(t)\ ) = 0. 



8\0 



s,t£B R+d(L) (0), \\s-t\\<5 



Representations of max-stable processes based on single extreme events 



27 



In particular, we can choose 5 > such that 

2|L© B R (0)\ " 



E sup |F( S1 ) - F(s 2 )\ < n . r Jl , n .. . (56) 



K s 1 ,s 2 &B R+d ( L) (0), 1 1 Si — s 2 ||<<5 

Then, we get 

pj sup \M n ( Sl ) - M n (s 2 )\ > e\ 

I — s 2 ||<<5, si,s 2 £i J 

<nPl sup |M n (si)-M n (s 2 )| >cne\ 

[ ||si— s 2 ||<(5, si,s 2 6l J 

<n(p{n(uu,t,f) : teL®B R (0), 

sup |/( S i)-/( S2 )|>— )) >0 

8i,8,eBR +<l (i)(0), u J / 

|]si-s 2 ]]<(5 

+ P : i GM d \(L©^(0)),sup|/(s-t)| > ^V) >0 

< n I 1 — exp < — / / 

V " I JL®B R (0) Je 

(CTIS \ 
sup |F(si) -F(s 2 )| > )ducdt 

I Sl — S2 \\<S 



+ 1 — exp 



d \(L®B«(0)) 



/ u _2 P ( sup|F(s - t)\ > — J ducdi 1 ) 
h(o))A Vsei u / J/ 



<nll — exp I E< sup — F(s2j 

V V n£ I Si,s 2 €B H+d(I) (0), 

s 2 ||<<5 

+ l-exp( / Eisup|F(s-t)|ldt 

V «e Jr\(l®b h (o)) UeL J 

< n (l- exp (-£)+! -exp (-£)) < „, 

where we used Equation (56) and (55). Thus, the sequence of processes {M n |x, : 
n e N} is tight. 
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